Analysis of the second order exchange self energy of a dense electron gas 
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Abstract 

We investigate the evaluation of the six-fold integral representation for the sec- 
ond order exchange contribution to the self energy of a dense three dimensional 
gas on the Fermi surface. 
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Introduction 



The second order exchange energy, represented by the diagram in Fig. la 
contributes importantly to the correlation energy of a dense electron gas [1] . It 
is given by the nine-fold integral 



E2x = 



327r4 



J dPpr j d^P2 J 



fpi fp2 fpi+qfp2+g 



^ {q + pi+P2)'^iq'^+pi-?+P2-q) 



(1) 



in three dimensions, where fp denotes the Fermi distribution function for elec- 
trons of wave vector p and denotes that for holes. In a remarkable display 
of mathematical virtuosity (1) was evaluated in closed form by Onsager[2] and 
Onsager, Mittag and Stephen [3] who found 



^2x = Jln(2)-^c(3). 



(2) 



Subsequently, Ishihara and Ioratti[4] worked out the corresponding value for a 
two-dimensional system, and the d-dimensional case was evaluated by Glasscr[5] . 

Recently the the second order exchange term in the electron self energy, 
represented by the diagram in Fig. lb was studied by Ziesche[6]. It is given, in 
three dimensions, by the six- fold integral 



q^ 



fpfk+qfp+ q fp fp 



pJp+q 



{k+p + q)^{(p + k- q+p- q) 



(3) 



For k = kp{= 1) Ziesche succeeded in decomposing (3) into the sum T,2x = 
— (Xi + X2)/47r^ of the two simpler integrals 



X2 = 



d^qi f d'^q2 fk+qi+q2fk+qifk+q2 



^2 



fqi 



qi ■ q2 

d^q2 /fe+gi+52'^'=+9i'^'=+92 

11 ^ 



(4) 

qi ■ q2 

and by following the procedure in [3], he managed to perform three of the 
integrations, thereby obtaining 



Xi 



-16n 



dp / dq 



dx 



F[p, q,: 



X2 = IStt / dp [ dq I 
Jo Jo J-l 



i(1-pV) l + q^ 

dx q'^Flp,q,x] 

(i-pV) 1 + 



(5) 



where 



a = 



1-q' 



2q 



2p 



a 



1 + p^q 



,2„2 



2pq 



F\p,q,x] 



■ tan 



ax + P 



(6) 
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The integrals in (6) are suitable for numerical evaluation and Ziesche found 
Xi = -30.70598 . . .,X2 = 21.28490 . . .. 

According to the Hugenholtz-van Hove- Luttinger-Ward theorem[7] T,2x = 

E2x, SO ^ 

Xi+X2 = 3C(3)-^ln(2). (7) 

The aim of this note is to evaluate X = X2 — Xi analytically, so as to obtain 
closed form expressions for the integrals in (4). 



Calculation 



From (5) we have 



X = 167r dq t dp t ^a: f^''l'^J . 
Jo Jo J-i 



(8) 



Since the limits on the x-integral are symmetric, we retain only the even part 
of the integrand of (8) by averaging X and the integral obtained by a; — > —x 
and combining the two arctangents, thus obtaining 



X = IGtt / dp dq 



tan 



dx- 



q2_/32 + i_j.2 



(9) 



Next, we set g = e ^, p 
a = cosh(u + v), and 

f.TT/2 



(1 — p'^q'^){a? — x"^) 
e~", X = sint/), so a = sinh u, P = sinh v 
X = 



(sinh(u+v)+sinh(ri— w)) cos 
sinh(tt+i;) sinh(u— v)+cos^ 



/•OO /'OC /> 

Stt / du dv d(p COS ^ 7 rf 27 ^ o — ^ 

Jo Jo Jo sinh(u + v) [sinh {u + v) + cos^ <p\ 



(10) 

/o ^0 

We make the coordinate transformation r = v+u, s = v—u, having Jacobian 
1/2, to obtain 



X = 47r 
Since 



dr / ds 



tt/2 tan 
dij) COS0 



(sinh r+sinh s) cos (p 
cos"^ (f>— sinh r sinh 5 



tan 



sinh r(sinh^ r + cos^ (/)) 
cos (^(sinh r + sinh s) 



(11) 



cos^ 



sinh r sinh s 



Im ln[(cos (f) + ismh r)(cos 0-|-isinh s)] 



tan 



_i f sinh 
\ cos ( 



(11) becomes 



+ tan 



X = 



_i f sinh s 
\ cos (j) 



(12) 
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47r / rfr / ds dcj) 

Jo J-r Jo 



COS 



^tan -"^(sec (/)sinh r) + tan ^(scc (f>smh s) 
sinh r(cos2 (j) + sinh^ r) 



(13) 

Once again, we may drop the term in the integrand of (13) odd in s and perform 
the elementary s— integration, so that 



X = 8tt 



rdr /■■^/^ _i / sinh r 

dcp tan ' 



sinh r Jq 



\ cos 



COS cf) 



cos^ (b + sinh r 



(14) 



To evaluate the (/)— integral, wc set tan ijj = sec ^sinh r , ^ = tan -"^(sinh r) 
cos~^(secft, r), to transform (14) into 



X = 87r 



rdr 



■ cos /i 



'^/^ cos tj) dij) 



sj sir? ijj ■ 



sin^ fj, 



sinh r 

The tjj— integral is tabulated [8] and X is reduced to a single integral 

2 f°° r sech r ln(l + sech r) 



X = 4r 



J 

Jo 



sinh r 

To evaluate the remaining integral, let 

r ln(l — a sech r) 
Q sinh rcosh r 



-dr. 



fia) 



dr 



(15) 



(16) 



(17) 



for which /(I) = X/Att"^ and /(O) = 0. By differentiation with respect to a and 
partial fraction decomposition, we obtain 



f°° rdr _ 2^ /-^ rdr _ 1 T 
Jo sinh r Jo sinh 2r a Jo 



rsinh r 



1 



1 



cosh r cosh r + a 



• (18) 



The first two integrals on the right hand side of (18) are tabulated[9] and, after 
an integration by parts, we find 

df TT^ 1 Z""" 

{1- a^J-j- = —{2- a) - - a sech r)dr (19) 

da, 8 o, Jo 

The substitution u = sech r leads to another tabulated integral[10], giving 



da 



8a V 1 + a 



1 (cos ^ a)^ 
2a l-a2 



+ 



(20) 



which, with the substitution a = cos 6 yields 



X = 4Tr^ r ^da = nHn(2) + Att^ ^^[0^ - —{1 - cos(29))]. (21) 
Jo da ^ ^ Jo sin 26*^ 8 ^ ^ ^ ^ 
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Finally, we find by setting (j) = 26, and folding the new range of integration 
[n/2,TT] back to [0,7r/2] 



X = 7r^ln(2) + 47r2 / 
Jo 



A(t>{(t> ~ it) 



sm 



7r4ln(2)--7r2C(3), 



where we have used[ll] 

^^/^ <t>d<t> 
iQ sin (f) 



f 

Jo 



= 2G, 



L 



sin (j) 2^^ ' 



(22) 



(23) 



in which G denotes Catalan's constant. 



Discussion 

Our result is that we have obtained closed form expressions for the two 
six-fold integrals in (4) 



ln(2)-^C(3) 



(24) 



-30.70598523924889925762268444608481536875855208165945918981645846 . . . 



Xo 



^ln(2)-Ac(3) 



(25) 



21.284905670516337983402598547497784400625730440810132220995696061 . . . 
This gives the value 

S2x = (26) 

0.0241791589181444058954507621628984314049152384251207335945309986 . . . 

in agreement with Zieschc's [6] seven place calculation. We hope to extend our 
calculation to an electron gas of arbitrary dimension, as was done for E2x- 



Acknowlegements 

The first author thanks Dr. Paul Ziesche for a discussion of his work and 
the National Science foundation for support under Grant DMR 0121146. 



5 



References 



[I] M. GcU-Mann and K. Brueckner, Phys. Rev. 106, 364 (1957). 
[2] L.Onsager, Unpublished. 

[3] L. Onsager, L. Mittag and M.J. Stephen, Ann. Physik (Leipzig) 18, 71 

(1966). 

[4] A. Isihara and L. loriatti, Phys. Rev. B22, 214 (1980). 
[5] M.L. Glasser, J. Comp. Appl.Math.lO, 293 (1984). 

[6] P. Zicschc, ArXiv.org cond-mat/0605188. (Annalcn der Physik, In Press.) 

[7] J.M. Luttinger and J.C. Ward, Phys. Rev.118, 1417 (I960). 

[8] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series and Products, 

6th ed., Academic Prcss(2000). p. 466, No. 3.842(2). 
[9] Ref.8, p.369, No.3.521(l) 
[10] Ref. 8, p.554, No.4.292(5). 

[II] Ref. 8, p. 427, Nos. 3.747(1,2). 



6 



